Transverse instability of dunes 
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The simplest type of dune is the transverse one, which propagates with invariant profile orthogo- 
nally to a fixed wind direction. Here we show numerically and with a linear stability analysis that 
transverse dunes are unstable with respect to along-axis perturbations in their profile and decay on 
the bedrock into barchan dunes. Any forcing modulation amplifies exponentially with growth rate 
determined by the dune turnover time. We estimate the distance covered by a transverse dune before 
fully decaying into barchans and identify the patterns produced by different types of perturbation. 
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Wind directionality and sand availability are the main 
factors dictating dune morphology. Bimodal and mul- 
tidirectional wind systems form longitudinal and star 
dunes, respectively [1]. Under unimodal winds, two types 
of dune may occur, depending on the amount of sand: 
crescent-shaped barchans, evolving on the bedrock, and 
transverse dunes, which appear when the ground is cov- 
ered with sand [l|43|. Studies of dune genesis have fo- 
cused on the growth of longitudinal sand- wave instabili- 
ties of a plane leading to transverse dunes — which mi- 
grate downwind with invariant profile orthogonal to the 
transport direction 0]. The stability of the transverse 
dune shape, however, has remained a long-standing open 
issue, of relevance for several areas of aeolian research 
and planetary sciences. As shown in water tank experi- 
ments, a transverse sand ridge of finite length evolving on 
the bedrock destabilizes and decays into barchans when 
subjected to a stream of nearly constant direction [5]. 

The complete quantitative study of transverse dune 
evolution requires a mathematical modeling that com- 
bines the description of the average turbulent wind field 
with a model for sand transport in three dimensions . 
Here we adapt this model in order to investigate system- 
atically, for the first time, the stability of a transverse 
dune under unidirectional wind. The dune model con- 
sists of iteratively performing the calculations listed in 
the steps which follow. 

(i) Wind — the average wind shear stress field (r) over 
the terrain is calculated from the equation, 

T= |To|(To/|To|+f), (1) 

where tq is the wind shear stress over the flat ground, 
and the shear stress perturbation due to the local topog- 
raphy, f, is computed by solving the three-dimensional 
analytical equations of Weng et al. [l^. Since this 
wind model is only valid for smooth surfaces, the cal- 
culation must be adapted in order to account for flow 
separation at the dune brink. For each longitudinal slice 
of the dune, a separation streamline, s{x,y)^ is intro- 
duced at the dune lee, where x and y are the direc- 



tions longitudinal and perpendicular to the wind, respec- 
tively. The wind model is solved, then, for the envelope 
hs{x^ y) = max{/i(x, y)} comprising the dune sur- 

face, h{x^y)^ and the separation streamlines at the dune 
lee; these define the so-called separation bubble, inside 
which the wind shear is set to zero [7]. The shape of 
s{x^y) is approximated by a third-order polynomial, the 
coefficients of which are calculated from the continuity of 
/i, s and their respective first derivatives at the brink and 
at the reattachment point downwind, which is computed 
assuming that s{x,y) has a maximum slope [7]. 

(ii) Sand flux — next, the mass flux of particles in salta- 
tion — which consists of grains travelling in ballistic tra- 
jectories and ejecting new particles upon colHsion with 
the bed [2] — is computed from a continuum model 0, 
in which the saltation cloud is regarded as a thin fluid- 
like layer moving over the immobile sand bed. When 
the wind shear stress exceeds a minimal threshold and 
saltation begins, the sand flux first grows exponentially 
due to the multiplicative process inherent to saltation 
transport. However, since the grains accelerate at cost 
of aeolian momentum, the flux cannot increase beyond a 
maximum value. This flux is reached after a saturation 
transient where the air shear stress within the saltation 
cloud equals the minimum value (rt) for sustained salta- 
tion [I, llH . This mechanism of flux saturation is taken 
explicitly into account in the calculation of the height- 
integrated mass flux per unit length and time (g), 

V-q = (l-|q|/9s)|g|/4, (2) 

where = [2a\vs\/g]{r — n) is the saturated flux, 4 = 
[2a\vs\'^ / g^]rt{T — Tt)~'^ is the characteristic length of 
flux saturation, Vs is the steady-state grain velocity, g 
is gravity, and a ^ 0.4 and 7 ^ 0.2 are empirically de- 
termined model parameters p. 

(iii) Surface evolution — the local height is updated from 
the equation, 

dh/dt = - V-gr/psand, (3) 
where Psand is the sand bulk density. Wherever the lo- 



2 



cal slope exceeds the static angle of repose of the sand 
34°), the surface relaxes instantaneously through 
avalanches in the direction of the steepest descent [7|. 
Eq. ([3|) is then iteratively solved, using the flux of 
avalanches along the slip-face, 

gavai = /c[tanh(V/i) - tanh(l9dyn)]V/i/|V/i|, (4) 

where A: = 0.9 and ^dyn = 33° is the dynamic angle of 
repose, until the local slope is below Odyn- 

The calculations are performed with constant upwind 
shear stress (tq), and open and periodic boundaries in the 
directions x and respectively. First, a transverse sand 
ridge of Gaussian cross section and invariant profile or- 
thogonally to the wind is let to evolve, under zero influx, 
into a transverse dune of height H and width Lq. The 
dune propagates, then, with fixed profile ho{x^y). The 
fragmentation of the dune into an array of barchans, re- 
ported from experiments [5|] , is not observed in the simu- 
lations. This result is independent of the amount of sand 
on the ground: the dune never breaks into barchans if 
Ty = dho/dy = 0, i.e. if the wind is unimodal and there 
is no variation in the transverse profile of the dune. 

Next we add to the dune profile a small perturbation 
of the form h{y) = ^o^(^), where (j){y) = cos[27r?//A], 
with 6o/H <C 1 and A constant, such that the modified 
dune profile reads h{x,y) = ho{x^y) + h{y). As shown 
in Fig. la, the perturbed dune is unstable, and decays 
after some time into a chain of barchan dunes, no mat- 
ter the values of Sq or A [12]. Differently from the be- 
havior of transverse instabilities of flat granular surfaces 
forced by an initial modulation [l3|, there is no thresh- 
old wavelength for the growth of the perturbations. The 
transverse dune is unstable also when (j){y) is a random 
function of y (see Fig. la), or when the modulation is on 
the X?/— plane. 

The magnitude of the perturbation, S{t)^ defined as 
the difference between the maximum and the mininum 
height measured along the dune crest, increases expo- 
nentially in time, S{t) ~ exp(<jt). The same behaviour 
is found for the difference between the maximum and 
minimum mass of longitudinal slices (Fig. lb). No de- 
pendence on the parameters of the initial modulation is 
observed for the characteristic time scale 1/a [HI , which 
is found to decrease with the wind speed and to increase 
with dune size. From dimensional analysis, a scaling of a 
with Qo/H^ is expected, where Qo = Qsi'To) / Psand is the 
saturated bulk sand flux associated with tq. The best fit 
to the simulation data leads to the expression (Fig. Ic), 

cT^0.32Qo/H^ = 1/T^, (5) 

i.e. the time scale of the instability growth is in essence 
the turnover time, Tm, of a barchan dune that has the 
same height H as the initial transverse dune Q. 

In this manner, dune genesis involves two different 
kinds of sand- wave instabilities. The longitudinal one. 




FIG. 1. (a) Spatio-temporal sketch showing the dune pro- 
file at different times plotted at the corresponding distance 
of migration downwind. Dune height is iif ~ 4.0 m and 
the wind blows from the top with constant shear stress 
To/rt ~ 3.1. The initial perturbation is random with ampli- 
tude ^0 = 12.5/xm, i.e. equal to the surface roughness (zo), ap- 
proximately 1/20 of a grain diameter 6]; the dune's turnover 
time is Tm ~ 6 x 10^ s. (b) Time evolution of the modulation 
amplitude, m(t) — 1 - Mmin/Mmax, where Mmin and Mmax 
are, respectively, the minimum and maximum mass of the 
longitudinal slices at time t; m{t) increases as exp(crt). (c) 
a scales with Qo/H^ (Eq. ([5|)). The continuous line, which 
has slope ~ 0.32, is the best fit to the simulation data us- 
ing Eq. (O. Circles: simulations with ro/rt ~ 2.6 and 3 m 
< if < 18 m; squares (triangles): H — 4.0 m (6.0 m) and 
1.5 < To/rt < 6.0. A simulated barchan, with central slice of 
length Lo and corresponding cross-wind width VK, is shown 
in the bottom right-hand corner. 



which leads to the transverse dune 0, 0, results from 
the combined effect of an upwind shift of the shear stress 
with respect to topography and the downwind lag of sedi- 
ment flux with respect to flow — the so-called saturation 
length (cx 4), which is proportional to ^drag, i-e. the 
grain diameter multiplied by the grain to fluid density 
ratio |l5[ . The transverse instability, studied for the 
first time in the present work, is consequence of cross- 
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FIG. 2. Unstable patterns of a transverse dune of height 
H = 3.0 m under To/n ~ 3.1. Snapshots of simulations ob- 
tained with A = 4 m (a), 64 m (b), 128 m (c) and random 
perturbation (d), all at \m\ ^ 99%. Each image has dimen- 
sions 512 m X 110 m. The plots show the corresponding 
auto-correlation functions C{y) = {m{y + y')m{y')) for the 
residual mass profiles m{y,t) = M(y) — {M(y)), where M(y) 
is the mass of the y—th longitudinal slice and C(0) is the 
square of the standard deviation. The first maximum of C{y) 
gives the average wavelength. The arrow indicates the width 
(W ~ 41 m) of a barchan that has the same height as the 
transverse dune. 



wind sand transport along the dune axis, which plays 
a major role in coupling the longitudinal slices of dunes 
[l,^. We find from a linear stability analysis (see Section 
2 of the Supplementary Material ^]) that, since differ- 
ent slices can have different velocities due to the relation 
V ^ the lateral transport becomes the destabilizing 
factor for tranverse dunes. 

In fact, when the perturbation is small and the lateral 
wind component is negligible, cross-wind transport oc- 
curs mainly due to gravitational downslope forces arising 
wherever on the slip- face the local slope exceeds the angle 
of repose [11,0] • It can be easily verified that mass transfer 
between neighbouring slices occurs if there is a relative 
shift in their downwind positions: sand moves from the 
slices advanced downwind toward the retarded ones . 
The relation v ^ 1/H implies that the smaller slices, 
which migrate faster, move advanced with respect to the 
larger ones. Thus, if a small perturbation is applied to 
the mass profile of a transverse dune, then, due to lateral 
transport, the smaller slices lose mass and become even 
smaller and faster, while the larger slices become even 
larg er and slower, such that the perturbation increases 
[12|. It can be shown that, in this linear regime (i.e. for 
small perturbation), the growth rate of the perturbation 
scales with 1/A [12]. As the perturbation increases, non- 
linear effects become important and the study of dune 
evolution must rely on numerical simulations (Fig. 1), 
which will be discussed next. If the dune is on the bed- 
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FIG. 3. Unstable dune patterns occuring in nature, to- 
gether with the corresponding patterns obtained in the cal- 
culations. (a,b) dunes at White Sands and (c) Guerrero 
Negro, Baja California (images courtesy of Labomar); (d) 
barchans emerging from an unstable transverse dune in a wa- 
ter tank experiment ^] (image credit: Sylvain Courrech du 
Pont - Laboratoire MSG, Univ. Paris Diderot, Paris, France); 
(e) Martian barchans near 7.4°N, 292. 3°W (image courtesy 
of NASA/JPL/MSSS). 



rock, then, when the perturbation becomes of the order 
of the dune height, barchans separate as sand is released 
through the emerging limbs along the barchanoidal chain. 

The pattern of the emerging barchan chain depends on 
A (Fig. 2). When A is approximately W 12H + S^drag 
[9], i.e. the width of a barchan that has the same height 
{H) as the transverse dune, a nearly uniform chain of 
barchans emerges, all of width W. For larger A the 
dune destabilizes in spatial cycles of wavelength A, which, 
again, decay into smaller barchans. Periodicity is lost 
when A < ly or when the perturbation is random. Never- 
theless, at the time when barchans are about to separate, 
the characteristic wavelength of the unstable pattern in 
all cases is of the order of as can be seen in Fig. 2 from 
the auto-correlation function of the mass profile in trans- 
verse direction. Interestingly, simulations using the afore- 
mentioned parameters can qualitatively mimic different 
dune patterns found in Nature, as depicted in Fig. 3, 
suggesting that the patterns are quite generic. Our sim- 
ulations show that all such patterns evolve toward an 
array of separated barchans, independently of the dune's 
past development stages or on the physical conditions. 

From the expression S{t) ~ Joexp(t/Tm), it follows that 
a transverse dune of height H and profile modulated by a 
perturbation along the axis of initial amplitude do needs a 
time Too ^ Tm^n{H/So) to fully decay into barchans. The 
total length Loo covered by the transverse dune during 
the decay process should scale, thus, with W\n{H/So). 
Using the relation W ^ cLq (c.f. inset of Fig. Ic), where 
c ^ 2.82 • [0.88-Tt/ro] for ro/n > 2 [12], we obtain 
Loo ^ cLoln{H/So). Since 6o cannot be smaller than the 
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roughness of the sand surface (zq), an upper bound for 
Loo can be estimated from the expression, 

Loo <cLoln(i7/zo). (6) 

By assuming zq ^ 10/im 0, [Tej and c ^ 1.5 0, i, S, [I^, 
a transverse dune with H ^ 4 m and Lq ~ 29 m should 
fully break into barchans within a distance shorter than 
Loo ^ 19Lo ~ 550 m. Since all information on wind 
speed and on the attributes of sediments and atmosphere 
are encoded in Tm and in the dune's morphological rela- 
tions, Eq. is universal, i.e. it can be used to predict 
Loo of a transverse dune under any physical condition. 

In fact, the prediction of Eq. (|6j) is in agreement with 
observations from water tank experiments on transverse 
dunes moving under constant water stream (see Suppl. 
Mat. [l2| for an image). In the experiments [5|, a sand 
bar subjected to a flow of nearly constant direction is 
shaped into a transverse dune of height H ^ 2 mm and 
width Lo ~ 3 cm. The dune destabilizes and gives rise 
to barchans, which separate after a distance Loo ^ 3Lo 
Q. Indeed, by taking H = 2 mm and the value c = 
W/Lq ^ 1 of the barchans produced in the experiments 
jsf, Eq. ([6]) predicts Loo ^ 5.3 Lq, which is close but 
above the experimental value. Since in real conditions the 
initial perturbation is normally larger than zq observed 
values of Loo should be, indeed, always smaller than the 
theoretical estimate (Eq. (|6])). Moreover, varying wind 
directions are also a destabilizing factor and can further 
accelerate the decay process of transverse dunes [HI, |l3 • 

In summary, our results show that transverse dunes are 
unstable with respect to along-axis perturbations in their 
profile and decay into barchans if moving on the bedrock. 
Any instability amplifies with rate 1/Tm, whereas the 
cross-wind width of the emerging barchans scales with 
the dune height. Our calculations show that a trans- 
verse dune emerging, for instance, from a sand beach in 
a coastal area should, after reaching the bedrock, migrate 
only a few times its own width until decaying into a chain 
of barchans — indeed, small barchans can be observed 
already when the first dunes leave the sand beach and 
enter the bedrock area where the ground is not covered 
with sand [l8|. These findings are clearly important to 
understand the mechanisms of dune size selection and 
the genesis and evolution of barchan dune fields. 

The transverse instability explains why barchans are 
the dominant dune shapes in areas of low sand avail- 
ability ^, i3i]. Our linear stability analysis \i2^ in fact 
shows that any transverse dune is always unstable. Since 
the linear stability analysis cannot go beyond small de- 
viations from the perfect translational invar iance, it can 
not discriminate between a wavy shaped transverse dune 
with eventually propagating waves, a barchanoidal chain 
or a decomposition into individual barchans. So various 
scenarios might be possible if the ground is not a bedrock 
but mobile, wet or vegetated. This open question might 
be tackled in the future with large scale computations. 



Transverse dunes are in fact also unstable under vari- 
ations in wind direction [HI, or collisions with other 
dunes [19]. While field observations are plagued with 
uncontrolled weather conditions and large time scales, it 
would be interesting to perform systematic laboratory 
experiments of the instability growth in order to confirm 
the results of our calculations. 
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